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ABSTRACT 

We explore the connection between the local escape velocity, V csc , and the stellar population 
properties in the ATLAS 3D survey, a complete, volume-limited sample of nearby early-type 
galaxies. We make use of ugriz photometry to construct Multi-Gaussian Expansion models 
of the surface brightnesses of our galaxies. We are able to fit the full range of surface bright- 
ness profiles found in our sample, and in addition we reproduce the results of state-of-the-art 
photometry in the literature with residuals of 0.04 mags. We utilise these photometric models 
and SAURON integral-field spectroscopy, combined with Jeans dynamical modelling, to de- 
termine the local V csc derived from the surface brightness. We find that the local V osc is tightly 
correlated with the Mg b and Fe5015 linestrengths and optical colours, and anti-correlated with 
the H/3 linestrength. In the case of the Mg b and Colour - V csc relations we find that the relation 
within individual galaxies follows the global relation between different galaxies. We intention- 
ally ignored any uncertain contribution due to dark matter since we are seeking an empirical 
description of stellar population gradients in early-type galaxies that is ideal for quantitative 
comparison with model predictions. We also make use of single stellar population (SSP) mod- 
elling to transform our linestrength index measurements into the SSP-equivalent parameters 
age (t), metallicity ([Z/H]) and a-enhancement [a/Fe]. The residuals from the relation are cor- 
related with age, [a/Fe], molecular gas mass and local environmental density. We identify a 
population of galaxies that occur only at low V CS c that exhibit negative gradients in the Mg b - 
and Colour - V GS c relations. These galaxies typically have young central stellar populations and 
contain significant amounts of molecular gas and dust. Combining these results with N-body 
simulations of binary mergers we use the Mg b- V csc relation to constrain the possible number 
of dry mergers experienced by the local early-type galaxy population - a typical massive ETG 
can have experienced only ~ 1.5 major mergers before becoming a significant outlier in the 
Mg b-V csc relation. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: abundances - galaxies: formation 
- galaxies: evolution. 



1 INTRODUCTION 

Early-type galaxies (ETGs including both ellipticals; E and lentic- 
ulars; SO) in the nearby universe represent an advanced stage of 
galaxy evolution. They likely result from the interplay of a wide 
rang e of physical processes . In the context of the ATLAS 3D sur- 
vey dCappellari et alj|2011al hereafter Paper I) we found that ETGs 
are well separated in terms of their stella r angular momentum 
dKrainovic et alj|201ll ; Emsellem ^TalfcOllL hereafter Paper II and 
Paper III respectively). Slow rotator ETGs can be explained as hav- 
ing accreted a major part of their mass due to hierarchical merg- 
ing dKhochfar et alj|201ll . hereafter P aper VIII), while fast rotators 
likely experienced cold gas accretion jDavis et al J201 ll i lSerra et all 
|2012| . herafter Paper X and Paper XIII) but then had the ir gas 
stripped by the cluster environment dCappellari eTai]|2011bl here- 
after Paper VII) or AGN and supernova feedback. Their present-day 
properties - their stellar populations, dynamics, gas content and 
environment - contain information about their entire evolutionary 
history. Through a detailed study of these properties we can hope 
to disentangle some of this history and learn about the key physi- 
cal mechanisms that have played a dominant role in shaping these 
objects. 

Early-type galaxies have a few advantages when it comes to 
studies of this nature. They are massive and bright, allowing high 
quality data on both their stellar content and dynamical state to be 
obtained with relative eas e. Their light is typi cally dominated by 
old, evolved stars (see e.g. iThomas et alj|2005) - the low levels of 
star formation allow us to make the assumption that their stars were 
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formed in a single event. The morphology of their central regions 
is relatively featureless, allowing us to study how properties vary 
within these objects without the complication of substructure such 
as spiral arms or clumps of star formation. 

Galactic archaeology is the study of the 'fossil' record within 
local early-type galaxies. A detailed analysis of the stellar popula- 
tions of these objects can reveal much about when the stars that now 
comprise these objects formed, but it reveals little of where they 
formed. In contrast, dynamical studies can reveal some of the where 
(in situ vs. accreted) but little of the when. By finding clear links be- 
tween dynamical and stellar population quantities we can begin to 
answer all of these questions, forming a picture of when and under 
what conditions the stellar populations of early-type galaxies were 
formed and assembled. 

There are several well-studied relations linking stellar pop- 
ulation and dynamical quanti ties; the colour-magnitude relation 
dVisvanathan & Sandagelll977h and the Mg-<r relation. Most early- 
type galaxies lie on the tight 'red sequence' of the colour-magnitude 
diagram (e.g. lBower et al , 1998 ), though some lie in the 'green val- 
ley' transition region or even in the 'blue cloud' jS trateva et al. 
1200 ll ; IConselicd[200rj ; Ivan den Bergh]|2007l ; iBernardi et aT]|2010l) . 



The Mg-cr relation has been studied b y man y authors (ex amples 



include iBender et alj 1 19931; Ijorgensenl ll997l; IColless et alj 
iTrager et al .1 l2000l : IBernardi et alj 1200 3l : iKuntschner et alj 



1999: 
2006) 



with many hundreds of early-type galaxies and though the precise 
zero points and slopes vary between studies the small scatter com- 
pared to the range in velocity dispersion covered is universal. 

The two above relations describe a connection between differ- 
ent galaxies, however we can also stu dy the variation within indi - 
vidual galaxies. This was first done bv lFranx & Illingworthl d 19901) 
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Figure 1. Comparison between total r-band magnitudes derived from INT 
and SDSS DR8 images for the 18 ATLAS 3D galaxies included in the latest 
SDSS data release. The solid line shows a fit to the data and the dashed line 
shows a 1 : 1 relation. The data are consistent with this 1 : 1 relation, with a 
scatter of 0.17 magnitudes. 



who studied the correlation of the local colour with the local a 
and escape velocity, V csc . They found that the global relation be- 
tween different galaxies is the same as the l ocal relation within 
indi vidual galaxies. This was f ollowed up by iDavies et alj i 19931) 
and ICarollo & Danzigeil dl994h . who studied the connection be- 
tween the Mg line strength and the local V esc . These works were 
all based on small samples of galaxies and relied upon long slit 
data, both in the measurement of line-strengths and as the input 
for the dynamical modelling necessary t o derive V e3 c. Two di - 
mensional spectroscopy was first used by lEmsellem et alj jl996l) , 
who used IF U data to study t he Mg 6-V eS c relation in a single 
galaxy, M104. IScott et alj J2009l . hereafter S09) investigated the line 
stre ngth - V eB c relations f or the SAURON early-type galaxy sam- 
ple dde Zeeuw et alj|2002h . significantly increasing the sample size 
(to 48 objects) compared to previous studies and also fully utilis- 
ing IFU spectroscopy, which allows us to study how kinematic and 
stellar population properties vary in two dimensions. 

The SAURON survey was biased towards high mass galaxies 
by the survey selection criteria, which uniformly sampled galaxies 
in absolute magnitude. In this paper we extend our study of the role 
of the local V csc to the complete, volume-limited ATLAS 3D survey 
of early-type galaxies (Paper I). In Section[2]we present the sample 
selection and observations. In Section[3]we present our method of 
determining the local V csc from the photometric and spectroscopic 
data. In Section [4] we present the results of our analysis, showing 
the correlation of various properties of early-type galaxies with the 
local Vesc - In the same section we extend this analysis to the study 
of the stellar populations of our sample. Finally, in Section [5] we 
discuss our results in the context of cosmological and N-body sim- 
ulations of early-type galaxy formation and state our conclusions in 
Section[6] 
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Figure 2. Upper panel: g-i colour map of NGC4753. Dust filaments show 
prominently in the upper left quadrant of the figure. Lower panel: g-i vs. 
semi-major axis distance for individual pixels form the above colour map. 
The red line shows a fit to the underlying colour profile. Pixels lying signif- 
icantly above this line are corrected for dust extinction as described in the 
text. At large radii the colour of individual pixels is dominated by noise and 
any uncertainty in the sky subtraction. 



2 SAMPLE AND DATA 

The ATLAS 3D project (Paper I) is a volume limited complete sur- 
vey of all nearby early-type galaxies (those not showing spiral 
structure in the visible photometry) brighter than Mk = —21.5 
(M ~ 6 x 10 9 M a ) and within 42 Mpc. The sample consists of 
260 galaxies spanning a broad range in central velocity dispersion, 
CT e , environment and absolute magnitude. In this work we consider 
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Table 1. WFC INT filters used for ATLAS 3D photometry 



Serial Number Name Central Wavelength FWHM 

(A) (A) 



204 


WFCRGOU 


3581 


638 


220 


WFCSloanG 


4846 


1285 


214 


WFCSloanR 


6240 


1347 


215 


WFCSloanI 


7743 


1519 


195 


WFCRGOZ 


8763 





a sample of 256 objects, excluding 4 galaxies for which poor quality 
or missing data (see below for details) prevented a full analysis. 

All 260 galaxies in the sample were observed with the 
SAURON integral-field unit (IFU) on the William Herschel Tele- 
scope at the Roque de los Muchachos observatory in La Palma. 
The SAURON data consists of integral-field spectroscopy over the 
wavelength range ~ 4800 — 5300 A out to typically 1 effective ra- 
dius (R e ). This wavelength range includes the three Lick stellar ab- 
sorption line strength indices; Mgb, Fe5015 and H/3. Maps of the 
stellar kinematics (parameterized by the mean velocity v and ve- 
locity dispersion a) and line strengths were produced as described 
in Paper I and McDermid et al., (in preparation) respectively. The 
velocity maps are presented in Paper II. 

2.1 Photometry 

We have obtained ugriz photometry for 258 galaxies from the 
ATL AS 3D sample, taken fro m the Sloan Digital Sky Survey (SDSS) 
DR7 dAbazaiian et"alll2009l) where available, supplemented by our 
own observations with the Wide Field Camera (WFC) on the 2.5m 
Isaac Newton Telescope (INT), again at the Roque de los Mucha- 
chos observatory. Observing conditions prevented the acquisition 
of photometry for the remaining two galaxies. We also note that 
two other galaxies suffered from low signal-to-noise ratio in their 
SAURON observations (see McDermid et al, in preparation, for de- 
tails) and we exclude these. 

Observations with the INT WFC were carried out to obtain 
ugri-band imaging for galaxies not covered by SDSS DR7. The ob- 
servations were obtained in 3 runs: 7-9 May 2007, 3 & 4 November 
2007 and 13-17 May 2008. Part of the observations were obtained 
under non-photometric conditions. Images were taken through the 
5 filters for 55 galaxies from the ATLAS 3D sample. The filters that 
were used are listed in Table [T] Integration times were typically 
60 to 160 seconds, reaching sensitivities comparable to or deeper 
than the SDSS. Fifteen galaxies were already observed by SDSS 
DR5 and were used to perform cross-checks in general and to bring 
the INT imaging onto the same photometric system as SDSS in 
particular. A further eight galaxies were made available in SDSS 
DR7 - these were not used for the photometric calibration but re- 
placed the corresponding INT observations. The i mages were re- 
duced and calibrated using the Astro- WISE system l lValentiin et alj 
|2007|) . iMcFarland et alj d201ll) give a detailed description of the 
Astro- WISE pipeline. 

The observations were debiased and flatfielded using nightly 
biases and twilight flatfields. Hot and cold pixels are identi- 
fied and indicated in the weight maps of the science images. 
Photometric calibration was derived in the Sloan ugriz sys- 
tem f rom nightly standard star field observations using SDSS 
DR5 dAdelman-McCarthv et al.1 12007) photometry in those areas 
where the INT imaging overlapped the SDSS survey coverage. 




10 20 30 



Figure 3. Photometry of NGC2685, the Helix Galaxy, with MGE models 
overplotted before (upper panel) and after (lower panel) dust correction, as 
described in the text. The effect of dust is visible as the band-like structure 
in the upper left of the galaxy; this structure is absent in the dust-corrected 
image. 



For the g- and /-bands w e used a linear colour transformation 
dVerdoes Kleijn et alj[2007h to convert from the WFC instrumen- 
tal photometric system to the Sloan standard photometric system 
. They found no evidence for a transformation in other bands. We 
verified the photometric calibration by using galaxies with both INT 
and SDSS DR7 coverage by comparing magnitudes of stars found 
in both images. Astrometric calibration was p erformed using th e 
USNO catalog as astrometric refere nce system {Monet et al .120031) . 

In the SDSS DR8 data release jAihara etal.l201ll) imaging for 
18 further ATLAS 3D galaxies was made available. We used this 
new SDSS photometry to assess the absolute photometric calibra- 
tion of our INT photometry. We derived total r-band magnitudes 
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Figure 4. Example MGE models of barred galaxies. Contours of the best-fitting MGE models (red) for NGC3941, NGC4267 and NGC4643 plotted over the 
photometric contours (black). The region shown was selected to approximately cover out to double the extent of the bar. The MGE matches the contours both 
inside and outside the region where the bar dominates. These models were produced automatically using the method described in the text to avoid fitting the 
bar. For galaxies with very strong bars (such as NGC4643) the fitting method is unable to match the photometry in detail, however the overall shape and surface 
brightness of the object are well reproduced. 



from t he two sets of images using the Multi Gaussian Expansion 
(MGE lEmsellem et aljfl994T) models described in Section \$A\ The 
comparison is shown in Fig. [TJ The data are consistent with a 1 : 1 
relationship between the INT and SDSS derived magnitudes, with a 
scatter of 0.17 magnitudes, indicating our INT photometry is well- 
calibrated. 



3 MASS MODELLING 

In order to determine the local V csc we must first construct a dy- 
namical model of each of our galaxies from our photometric and 
kinematic data. This involves several steps, summarised here and 
where necessary described in more detail below. 

i) The r-band image is corrected for internal extinction due to dust. 

ii) Nearby stars and galaxies are masked, as are regions where the 
dust correction failed. 

iii) A calibrated Multi Gaussian Expansion (MGE) model of the 
photometry is constructed from the mas ked and extinction correc ted 
image using the method and software o f lCappellari et all d2002l). 

iv) A series of Jeans Anisotropic MGE (JAM) models dCappellaril 
120081) are constructed using the photometric model with different 
values of inclination, i, mass to light ratio, M/L and anisotropy, j3. 
The models do not explicitly include a dark matter halo and the 
measured M/L represents the total one, including both luminous 
matter and any possible contribution from dark matter. 

v) The best fitting JAM model is selected by minimising the x 2 dif- 
ference between the predicted second moments of t he velocity field, 
V« 2 + a 2 and the SAURON observations (see ICappellari et alj 
l2012bl hereafter Paper XIX, for details of the fitting process). 

vi) This JAM model constrains the M/L normalisation and is used 
to transform the photometric model into a mass model. The JAM 
model also specifies the gravit ational potential. In practi se, we make 
use of the models presented in ICappellari etall J2012al) . 

vii) For each bin in the SAURON maps the local, line-of-sight, 
luminosity-weighted V csc is calculated from the gravitational po- 
tential of the best fitting JAM model. 



3.1 Dust correction 

The presence of dust in a galaxy can significantly affect the accu- 
racy of our MGE models, both through the dimming due to extinc- 
tion and by altering the observed morphology of the galaxy. We 
attempt to remove this effect by correcting for this extinction and 
fitting the MGE model to the underlying surface brightness profile. 
This correction was applied to 25 of our galaxies, approximately 10 
per cent of the sample. 

We use the technique of ICarollo et alj d 19971) (who assumed 
the dust is a screen in front of the stellar emission), further devel- 
oped by ICappellari etafl {2002 ) to correct for the extinction due to 
dust. We first construct a g-i colour map of each galaxy. This is 
converted to a colour profile by plotting the g-i colour against the 
logarithm of semi-major axis distance along ellipses fixed to the 
global PA and e (from Paper II), log m, of each pixel in the image. 
We use a robust linear fit to this data, which minimises the influ- 
ence of dusty regions, in order to determine the underlying colour 
gradient of the galaxy. We subtract this best fitting colour gradient 
from all pixels in the g-i colour map to determine the colour excess, 
E(g-0- Using a cut in E(g-i) (typically as a function of log m) we 
select pixels which have been significantly affected by dust extinc- 
ti on and correct for thi s using the standard galactic extinction law 
of lSchlegelet~al]dl998l) . The r-band extinction, A r , is related to the 
E(g-i) colour excess by the equation: A r = 1.15 E(gr — i). This pro- 
cess is illustrated in Fig. [2] where in the upper panel we show the 
g-i colour map for NGC 4753 (with the dust filaments prominently 
visible) and in the lower panel we show the robust linear fit to the 
colour profile. 

While the extinction correction determined by this method is 
not perfect, it improves the MGE fit to the photometry in the ma- 
jority of galaxies with significant amounts of dust. An example of 
this improvement is shown in Fig. [3] for NGC2685. The upper panel 
shows the r-band image and MGE model fit without dust correction, 
the lower panel shows the MGE model fit to the dust-corrected im- 
age. In the two most extreme cases of the galaxies whose photom- 
etry was significantly affected by dust, NGC4710 and NGC5866, 
this procedure failed to provide a reasonable correction for the r- 
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band photometry. For these galaxies we masked the dust-obscured 
areas from the MGE fitting procedure. 



3.2 Multi Gaussian Expansion modelling 

Using the MGE modelling technique of lEmsellem et al. 

I Gail), 

we 

constructed models of the surfac e brightness of eac h galaxy in the 
sample. We used the procedure o f lCappellaril l [2002l) as the basis for 
our technique but embedded it into our photometric pipeline to fa- 
cilitate dealing with a much larger sample of galaxies. This involved 
automating the process as much as possible. The position angle 
(PA) and centre of each galaxy was determined from the weighted 
second-moments of the surface brightness above a given level using 
the find_galaxy IDL routing. The level was chosen to reflect the 
global behaviour of the galaxy and to avoid the effects of bars and 
other non-axisymmetric features in the inner regions. Foreground 
stars and nearby galaxies were masked by hand where necessary. 

The MGE fits to the r-band photometry were performed by 
keeping the PA of the Gaussians constant in order to produce an ax- 
isymmetric MGE model to be used with the JAM modelling tech- 
nique described below. Each MGE model was convolved with a 
Gaussian point spread function with a — 1.0 (consistent to the 
typical seeing of the SDSS observations), before comparison to the 
observations and selection of the best-fitting MGE parameters. The 
resulti ng MGE model s are a ll corrected for galactic extinction fol- 
lowing ISchiegeletal] dl998l) . as given by the NASA/IPAC Extra- 
galactic Database (NED). The MGE models were converted to a 
surface density in solar units in the SDSS ma gnitude system us- 
ing a n r-band solar magnitude, M r ,0 = 4.64 dBlanton & Roweisl 
120071) . 



The MGE models described in this work are used as input for the 
JAM modelling method, which calculates a prediction of the line- 
of-sight second velocity moments (vf oa ) for a given set of mode l 
parameters and fits this to the observed V rmB dCappellarill2008l). 
We ma ke use of the simplest set of models from Cappellari et al.l 
(2012a). model (A). These are self-consistent axisymmetric JAM 
models, in which the dark matter is assumed to be proportional to 

the stellar mass. 

The potential <E> is calculated as in lEmsellem et al.l d 19941) and 
V csc is simply related to this by V csc = ^/2j<3?(i?, z)\. In or- 
der to deproject our 3-dimensional potential to compare with our 
SAURON index maps we assume that V csc is related to the indices 
by a power law relation of the form: Index oc V2 SC - With this as- 
sumption we can extract the luminosity-weighted average V eS c,p of 
the local V esc along the line-of-sight and produce V csc maps across 
the entire SAURON field (see S09 for details). We emphasise that 
the quantity we derive here and describe as the local V csc is es- 
sentially derived from the stellar surface brightness with a global 
scaling determined by the dynamical modelling. Our V csc is likely 
close to the true V csc in our galaxies, given that stars likely dom- 
inate the centre of ETGs (Paper XIX). However the quantity we 
compute does not intend to represent the true V csc , which is em- 
pirically uncertain due to the possible contribution of dark matter. 
Our aim is to find a purely empirical relation between quantities 
that can be robustly measured on real galaxies from photometry 
and integral-field spectroscopy. What matters is that our procedure 
can be replicated on simulated galaxies, to test whether they behave 
like real galaxies. When comparing our V esc to simulated galaxies 
one should determine the local V esc from the stellar particles alone. 



Bars and other non-axisymmetric features 

The MGEs were fitted in such a way as to avoid following bars or 
other non-axisymmetric features in order to reflect the underlying 
mass distribution of each galaxy. In our fully automated procedure 
this was achieved by constraining the maximum (q ma x) and mini- 
mum (q m in ) axial ratios of the individual Gaussians contributing to 
the fit. Successive fits were carried out, using a narrower range of 
axial ratios in each each subsequent fit until the mean absolute devi- 
ation of the model at a given q m in , q-max increased by > 10 per cent 
over the previous step. This simple prescription proved effective in 
reproducing the surface brightness of the majority of galaxies at all 
radii while avoiding the effects of bars and other non-axisymmetric 
features. Only in a small number of cases of very strong bars was 
this prescription ineffective. In these cases the MGE fit was carried 
out by hand. Examples of the MGE models of several barred galax- 
ies are shown in Fig. [4] A detailed analysis of the recovery of M/L , 
p z and i for barred objects is presented in lLablanche et al] d2012l . 
hereafter Paper XII), based on N-body simulations of axisymmetric 
and barred galaxies. There we find that the M/L recovered from the 
JAM modelling is typically accurate to within a few per cent, but 
for face-on objects can differ from the true value by up to 15 per 
cent. 



3.3 Jeans Anisotropic MGE (JAM) Modelling 

The dynamical model s used in this paper were presented in 
ICappellari et alj d2012al) and described in more detail in Paper XIX. 

Available as part of the MGE package of Cappellari (2002) 



3.4 Tests of MGE Photometry 

While MGE models of galaxy photometry have generally only been 
used as input to dynamical models, we would like to emphasise that 
these models provide a robust and flexible description of galaxy 
surface brightness profiles. This is true across a wide range of mor- 
phological types and is not based on any assumption about galaxy 
structure. In this section we compare photometric quantities (both 
total magnitudes and surface brightness profiles) derived from the 
MGE models described in Section [3] to studies in the literature us- 
ing standard techniques to derive their photometric quantities. Total 
magnitudes were derived from the MGE models using the following 
expression: 



E 27r/ >. 



2 / 



(1) 



where the sum is over the N component Gaussians and I'j, aj 
and q'j are the peak intensity, dispersion and axial ratio of the jth 
Gaussian compone nt. The primes denote quantities as measured on 
the sky plane (see ICappellari et al]|2002l for details). The surface 
brightness at a given position on the sky is given by: 



£(x',2/') = ^/j'exp 



2a? V q? 



(2) 



The only previous data set in the literature that contains pho- 
tometry for the entire ATLAS 3D s ample is the 2MASS Extended 
Source Catalog djarrett et alj|2000h . This allows us to test the ac- 
curacy and detect possible outliers over the full sample. In Fig. [5] 
we compare the absolute r-band magnitude derived from our MGE 
models (and distance-corrected using the distances in Paper I) to 
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Figure 5. 2MASS absolute S'-band magnitude vs. r-K colour derived from 
our r-band MGE models. The solid black line shows the best-fitting relation, 
with the dashed and dotted lines indicating the 1 and 2.6 a bounds. The blue 
points indicate > 2.6er outliers from the relation which were not included 
in the fit. 
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SDSS ModelMog (mag) 

Figure 6. Comparison between total r-band magnitudes derived from our 
MGE models and from the SDSS DR8 photometric object catalogue. Both 
magnitudes were derived from the same imaging data. The solid line shows 
the best-fitting relation, with the dashed and dotted lines indicating the 1 and 
2.6 a levels. Blue points were excluded from the fit. 



the absolute K-band magnitude from 2MASS. Overplotted is the 
best-fitting relation derived from an linear fit to the datapoints that 
minimises the absolute residuals after rejecting 2.6a outliers (en- 
closing 99% of values for a Gaussian distribution). We measure an 
observed rms scatter of 0. 1 1 magnitudes. This is a firm upper limit 
to the true error and intrinsic colour differences must be significant. 

The comparison with 2MASS magnitudes is between two 
datasets utilising different imaging data in very different bands. We 
also compared our derived apparent r-band magnitudes to appar- 
ent r-band Mo delMags taken from the SDSS DR8 photometric ob- 
ject catalogue dAihara et alj|201 it) . We restrict the comparison to 
galaxies which have no warning flags in their SDSS photometry, 
yielding a comparison of 214 galaxies. The comparison is shown 
in Fig. [6] For total magnitudes derived from the same images the 
scatter is significant, after rejecting 2.6 a outliers the rms scatter 
is 0.17 magnitudes. There is a systematic trend for the magnitude 
of the largest galaxies (R e > 30") to be underestimated by the 
SDSS ModelMags, with respect to our MGE based r-band magni- 
tudes. The offset is ~ 0.5 magnitudes. This trend was partially (but 
not completely) add r essed in the SPSS D R8 photometric pipeline 
dAihara et all 120 1 lh . iBlanton et al.1 d201 1) introduce an improved 
sky subtraction method that solves this issue, however it has only 
been applied to SDSS mosaic images and not the photometric cat- 
alogue. A detailed inspection of the SDSS DR8 photometric cat- 
alogue ModelMags for ATLAS 3D galaxies reveals several further 
problems. Some galaxies are significantly affected by confusion, 
where multiple PhotoObjects are identified for a single galaxy (e.g. 
NGC4710). Other ATLAS 3D galaxies are entirely missing from the 
photometric catalogue, despite lying within the survey area (e.g. 
NGC4486). Finally, a number of magnitudes in individual bands 
show significant errors (for example the (g-i) colour of -5.8 for 
NGC3384). These issues appear most significant for the ModelMag 
values, however the PetrosianMag values also exhibit similar prob- 
lems. 

To place the observed rms scatters of 0.11 and 0.17 magni- 
tudes (from the comparison to 2MASS if— band and SDSS r— band 
magnitudes respectively) in con text, we compare i t to the set of pho- 
tometric comparisons made bv lChen et all d20ld) . Among various 
comparisons they present from the literature, the most accurate is 
the one between their determination of total magnitudes gr and the 



one bv ljanz & LiskeJ d2008h . based on the same p-band SDSS im- 
ages and basically the same curve-of-growth method. They find an 
rms scatter in the two determinations of 0.09 mag (0.036 dex). This 
scatter is likely to represent the best of what can be achieved, when 
only minimal differences in the methods play a role. When compar- 
ing t heir qt values to thos e from the state-of-the-art ACSVCS sur- 
vey dFerrarese et al.|[20o3) in the same band, they find that the scat- 
ter increases to 0. 19 mag (0.076 dex), and the same a pplies for their 
compa rison against the VCC g-band magnitudes of iBinggeli et al.l 
dl984l) . Contrary to those comparison, we are using two very dif- 
ferent photometric bands (r and K s ), where intrinsic colour differ- 
ences must be significant. This is confirmed by the fact that the g — r 
colour-magnitud e relation, for 43 o bjects in common, between our 
r-band MGE and lChen et J] d201fj) g-band curve-of-growth deter- 
minations (their gA values) has an rms scatter of 0.062 dex, which 
would indicate an error of 1 1% in each of the two determinations. 
We conclude that our MGE photometry is as accurate as that of 
other state-of-the-art surveys. Moreover we only detect two possi- 
ble outliers, while there is no evidence for possible calibration is- 
sues for the entire sample. We therefore assign an error of 10% in 
our M r . 

Finally, to demonstrate that we are able to not only derive to- 
tal magnitudes from our MGE models but full surface brightness 
profiles we compare the major-axis surface brightness profiles de- 
ri ved from our MGE mod els to the major-axis surface photometry 
of iKormendv et al.l d2009l) . There are 23 galaxies in common be- 
tween their morphologically selected sample of Virgo cluster galax- 
ies and the ATLAS 3D sample. Surface brightness profiles for all 
23 galaxies are shown in Appendix B. We restrict the compari- 
son to the region over which our MGE models are reliable; 1,5" 
< r < 2max((jMGE), where max(fjMGE) is the a (in arcsec- 
onds) of the largest Gaussian in the MGE model for each galaxy. 
To account for the difference between the bands we shift our sur- 
face brightness profiles by a constant value such that the mean 
sur face brightness o v er the region of comparison matches that of 
the IKormendv et alj (|2009) profiles. The agreement between the 
profiles is excellent, with mean residuals (after applying the con- 
stant shift) of ±0.04 magnitudes. We accurately reproduce the 
surface brightness profiles for the full range of ellipticities in the 
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Figure 7. The Index- V csc relations for the whole sample. Blue line is a fit to R c circular aperture values for all galaxies. Red line is a fit to the same aperture 
but for only those galaxies with V csc r c > 400 kms -1 . With all three indices the relation is much tighter above V osc > 400 kms -1 (logV osc > 2.6). Below 
this value the scatter is significantly increased, largely due to a population of galaxies with unusually low Mg b values in their central regions. 



© 2012 RAS, MNRAS 000. [711191 



The ATLAS Project - XXI. V csc gradients 9 



Table 2. MGE parameters for the deconvolved r-band surface brightness 
for NGC4570, as an example of the ascii MGE parameter files available to 
download from our website. 



lop- V 

° 3 


log Oj 


n' 


(L0 pc ) 


(arcsec) 




45198.69 


0.318 


0.72 


9747.41 


1.142 


0.80 


3741.94 


2.743 


0.65 


1947.36 


3.796 


0.72 


935.15 


9.019 


0.54 


45.97 


13.037 


0.80 


495.96 


19.671 


0.20 


255.63 


36.549 


0.23 


37.80 


60.620 


0.28 


6.07 


60.620 


0.80 



Notes: Column (1): Total intensity of the j'th Gaussian component. Column 
(2): Dispersion of the j'th Gaussian component. Column (3): Axial ratio of 
the j'th Gaussian component. 

iKormendv e7ai]j2009l) sample. In summary, our photometric quan- 
tities, both total magnitudes and surface brightness profiles, derived 
from the MGE modelling are of excellent quality. The full set of 
calibrated MGE models will be made available from our websit^l 
at the conclusion of our project. An example MGE model is given 
in Table [2]- the columns (described in the table caption) are the 
same as those available from our website. The MGE models for the 
full sample (overplotted as contours on the r— band imaging) are 
shown in Appendix A. 



4 RESULTS 

Using the SAURON line-strength maps and the V eS c maps de- 
scribed above we constructed profiles by taking the average within 
elliptical annuli. The ellipticity chosen was a global ellipticity mea- 
sured using the weighted second moments of the surface brightness 
and was taken from Paper II. S09 found that this choice of elliptic- 
ity minimised the errors on each point within a galaxy's profile, and 
that their results did not depend on the choice of aperture (circular, 
fixed or varying e or major axis). We find a stronger trend of in- 
creased errors for non-optimal choices of aperture (likely due to the 
lower signal-to-noise (S/N) of the ATLAS 3D IFU data compared to 
the SAURON survey), however in general we confirm this lack of 
dependence on aperture. We also find that the relations presented 
in this work are robust against the choice of aperture. The error on 
each point for the line-strengths is the rms sum of the measurement 
errors and the rms scatter within each annulus. Foreground stars 
and bad bins due to low S/N were masked in the line-strength maps 
before extracting the profiles. 
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Figure 8. Each histogram shows the Index- V CS c gradients of the individual 
galaxies determined from a robust linear fit to each galaxy profile (black 
histogram). The red and blue histograms show the same data for galaxies 
with V csc p,c greater and less than 400 kms -1 respectively. The solid verti- 
cal line shows the global gradient determined from a linear fit to the central 
R e aperture values for the whole sample, with the dashed lines showing the 
3cr errors on this global gradient. For the Mg b gradient histogram, the peak 
of the distribution for individual galaxy gradients coincides exactly with the 
global gradient, though there is significant scatter in the individual galaxy 
gradients. For the Fe5015 and H/3 gradient histograms there is a significant 
offset between the peak of the distribution of individual galaxy gradients 
and the global gradient. 



4.1 The Index- Vosc relations 

The Index- V osc relations are shown in Fig. [7] with the profile for 
each galaxy shown as a black line. The blue line shows a fit to the 
central R c circular aperture values for the whole sample. The red 
line shows a fit only to the galaxies with V csc ,r > 400 kms - 
(log V csc > 2.6). All three indices show a tight correlation with 

2 http://purl.org/atlas3d 
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Vcsc, though the scatter increases at lower Vesc. This is particu- 
larly evident in the Mg b-V csc relation, where a clear break occurs 
at Vcsc ~ 400 kms -1 , though the effect is present with all three 
indices. This increase in scatter is primarily due to a population 
of galaxies having negative Mgfe- and Fe5015-V CS c gradients. In 
the H/3-Vcsc diagram these galaxies are seen to have steeply posi- 
tive gradients, compared to the typically flat profiles of the bulk of 
galaxies in the sample. 

The behaviour of the individual galaxy gradients can be better 
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seen in Fig.[8] The histogram shows the individual galaxy gradients 
determined from robust linear fits to each galaxy profile. The solid 
vertical line shows the global gradient, determined from a robust 
linear fit to the line-strengths and V csc as measured within a circu- 
lar R e aperture. In the case of Mg b the individual galaxy gradients 
are consistent with the global gradient, though there is significant 
scatter in the distribution about the global value. This is not true 
for Fe5015 and H/3, where the peak of the distribution of individ- 
ual galaxy gradients is significantly offset from the global gradient. 
When we divide the sample into high (V csc > 400 kms - ) and 
low V osc samples we find that the low V esc sample is biased to neg- 
ative gradients. In the high V csc sample we find only 2.7 per cent of 
galaxies have negative gradients, however this rises to 15.9 per cent 
in the low V csc sample. If we remove these negative gradient galax- 
ies from the sample we no longer observe an increase in scatter at 
low V csc . For the remainder of this section we exclude the negative 
gradient galaxies from our analysis, however they will be discussed 
further in Section [5~T1 The Index - V csc relations determined from 
galaxies with Vesc (Re) > 400 kms -1 (red lines in Fig. |7J are given 
below: 

log Mg b = (0.40 ±0.02) log V csc + (-0.53 ± 0.04) ; a = 0.04 
logFe5015 = (0.16 ± 0.02) log V csc + (0.23 ± 0.03); a = 0.03 
log H/3 = (-0.29 ± 0.04) log V csc + (1.02 ± 0.06); a = 0.06 

Here, a indicates the robust rms deviation orthogonal to the fitted 
linear relation. The Mg& and Fe5015 relations are ~ 3cr steeper 
than those found in S09, with improved formal uncertainties. The 
relation for H/3 is consistent with the S09 result within the un- 
certainties. This is because the ATLAS 3D sample contains a large 
number of lower mass ETGs, a population not well sampled by the 
SAURON survey. 



4.2 The Colour- Vcsc relations 

lFranx&IllingworthNl990l) originally studied the Colour- V esc re- 
lations, but recent work has focused exclusively on the correlation 
with line-strengths. With the ATLAS 3D sample we have complete 
ugriz coverage of the sample, allowing us to re-examine their orig- 
inal relation. Colour profiles were extracted in the same fashion as 
the line-strength and V CS c profiles; median colours were measured 
within elliptical annuli on the sky-subtracted, dust-corrected and 
masked SDSS and INT images. In Fig. [9] we show the g-r and u-r 
- Vcsc relations. We selected these colours as the g- and r- bands 
have the highest S/N and hence show the tightest relation, while 
the w-band is more sensitive to young populations. We find the 
same behaviour as in the Mg b - V OS c relation; a tight correlation 
with increasing scatter appearing below V CS c ~ 400kms~ 1 . This 
increase in scatter is more evident with the u-r colour than with 
g-r. We also find the same local and global behaviour described 
above (see lower panel of Fig. |9), in that individual galaxy gradi- 
ents are consistent with the global relation determined from cen- 
tral values. Fitting to the R e aperture values for those galaxies with 
Vesc, Re > 400 kms - we find that the colours are given by: 

g - r = (0.26 ± 0.03) log V CS c + (0.08 ± 0.05); a = 0.04 (3) 
u-r= (0.80 ±0.08) log Vcsc + (0.47 ±0.14); a = 0.12 (4) 

Below Vcsc ~ 400 kms -1 we again find a significant increase in 
the fraction of negative gradient galaxies; 3.9 per cent at high V OS c 
increasing to 22.0 per cent below ~ 400 kms -1 . 



Table 3. Principle components analysis 







Age' 


[Z/H]' 


[a/Fe]' 


Eigen- 


% of 












value 


variance 


PCI 


0.680 


0.187 


-0.137 


-0.696 


1.87 


47 


PC2 


0.350 


-0.609 


-0.644 


0.304 


1.38 


34 


PC3 


0.522 


0.554 


0.065 


0.646 


0.59 


15 


PC4 


0.379 


-0.537 


0.750 


0.078 


0.16 


4 



Notes: The primed variables are standardised versions of the corresponding 
variables with zero mean and unit variance. The coefficients of the principal 
components are scaled to the variance and sensitive to the range of each 
variable, in the sense that variables that only vary by a small amount tend to 
have a large coefficient. 



4.3 Vcsc and the globally averaged SSP parameters 

Line strengths and colours can be indicative of trends in the stellar 
populations of our target galaxies but the information we are really 
interested in is the star formation histories of these galaxies. These 
histories will typically be quite complex, involving multiple peri- 
ods of star formation over the course of a galaxy's lifetime with 
different durations, contributions to the total stellar mass and chem- 
ical compositions. As an initial test we can simplify this scenario 
by treating all our galaxies as single stellar populations (SSPs), i.e. 
all the stars were formed in a single, instantaneous burst, and we 
can characterise these SSPs by three parameters, the age, metal- 
licity and enhancement of alpha elements. While the SSP scenario 
is a simple one it can provide a useful starting point for studying 
early-type galaxies, where a significant fr action of their stars w ere 
typically formed at high redshift (see e.g. iThomas et alj|2005l) . In 
these cases the approximation of a single star-formation episode is 
a reasonable one. In particular it is important not to place too much 
emphasis on the precise values of the SSP parameters of individual 
objects, but the SSP models can be used to study the differences 
between objects and any trends arising from this. 

W e adopt the single stellar population models of ISchiavonl 
d2007l) to interpret our global line-strength measurements of Mg b, 
Fe5015 and H/3 in terms of the SSP parameters, age (t), metallicity 
([Z/H]) and alpha enha ncement ([a/Fe]). We fo llow the approach 
previo usly descr i bed by iKuntschner et alj ( l2010t ). where the mod- 
els of ISchiavonl 420071) were found to provide physically realis- 
tic SSP-equivalent ages, metallicities and alpha enhancements for 
the full range of line strength index m easurements found in that 
work. In particular, the lSchiavonl J2007h models were found not to 
satura te in age at high values of the H/3 index. IKuntschner et alj 
conducted a de tailed comparison between the predictions 
of the ISchiavonl d2007h models and other competing SSP models, 
finding their derived parameters did not depend significantly on the 
choice of SSP model. Further details of our global stellar population 
modelling are given in McDermid et al. (in preparation). While the 
determination of each of the three parameters depends on all three 
indices, some parameters depend more heavily on one index over 
another. H/3 can be broadly thought of as an age indicator whereas 
Fe5015 and Mg b are more strongly influenced by [Z/H] and [a/Fe]. 

Rather than examining the individual SSP-V esc relations (for 
a discussion of the SSP-V CS c relations see S09) we consider the 
four-dimensional parameter space defined by V C sc, age, [Z/H] and 
[a/Fe]. To analyse these data we use the technique of Principal 
Component Analysis (PCA). A description of the P CA technique 
and the interpretation of i ts output can be found in iFabed i 19731) 
and lFrancis & Wills! d 19991) . 
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Figure 9. Upper panels: The u-r and g-r - V csc relations. As before individual galaxy profiles are shown as the black solid lines with the global relation, 
again determined from R e aperture values, plotted in red. Both relations show the same local and global behaviour as the Mg fe-V csc relation. We again see an 
increase in scatter below V osc ~ 400kms due to a population of negative gradient, 'blue core' galaxies. Lower panels: Histogram of the individual galaxy 
gradients for the g-r - and u-r - V csc relation, determined by robust linear fits to each galaxy profile. The global relation is shown as the solid vertical line with 
the 3 cr errors shown as vertical dashed lines. The peak of the distribution of individual galaxy gradients coincides with the global gradient, and the width of 
the distribution is consistent with the error on the global gradient. The n-r-V osc relation also exhibits the same local and global connection. 



Applying PC A to our data we find that ~ 81 per cent of 
the scatter is confined to the first two principal components (Table 
[3), indicating that in this four-dimensional space local early-type 
galaxies are effectively confined to a plane. While this is consis- 
tent with the findings in S09 and the weights associated with each 
physical variable in the corresponding PCs are very similar, we do 
find a significantly increased scatter associated with the less signif- 
icant third and fourth principal components - we find a 'thicker' 
plane than in S09. Part of this effect can likely be attributed to the 
somewhat poorer quality of the IFU data in the ATLAS 3D survey 
compared to the SAURON survey due to i) 1 hour instead of 2 hour 
exposures and ii) a higher proportion of faint galaxies. Large values 
of PC3 and PC4 are typically associated with extreme values of age 
or [a/Fe], lying on the boundaries of our SSP model grid. Interest- 
ingly, the negative gradient galaxies do not occupy any special part 
of the parameter space and, in this respect, are indistinguishable 
from typical early-types. 

We take this idea further by attempting to fit a single plane to 
the local V osc , age and [Z/H] data points for our full sample. Using 



a robust, linear plane fit to our data we find that local early-type 
galaxies can be described by the relation: 



log 



V e 



V 500kms~ 



= 0.62 



-0.26 log 



( 1 



VGyrs 



-0.17 (5) 



with a = 0.07. This relation is shown in Fig.[l0] This SSP - V csc 
relation is much tighter than any of the individual age, [Z/H] or 
[a/Fe] - Vcsc relations. However, this relation shows an offset be- 
tween the global gradient (0.63 ± 0.03) and the distribution of the 
local gradients (0.75 ± 0.51, see the lower panel of Fig. Hot. If we 
use the relation previously found by S09 we find a somewhat in- 
creased scatter of a — 0.11 (though still significantly tighter than 
the individual SSP-V OS c relations) but we find a better match be- 
tween the local (0.97 ± 0.04 ) and global (0.95 ± 0.71) gradients. 
We also find that this SSP - V csc relation displays the same local 
and global behaviour as the Mg b- and Colour- Vesc relations. 
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S3P-V eso gradient 

Figure 10. Upper panel: The best-fitting SSP - V csc relation determined by 
robustly fitting a plane to the local V esc , age and [Z/H] of our galaxies. As 
before, the solid black lines show individual galaxy profiles and the red line 
is a robust linear fit to the R e aperture values for the whole sample. Lower 
panel: The distribution of the individual galaxy gradients. The solid vertical 
line shows the global gradient determined from the R e apertures and the 
dotted lines show the 3-sigma errors on this gradient. 



4.4 Deviations from the Mg b-V OS c relation 

Because of the tightness as well as the global and local behaviour of 
the Mg b-V C s C relation we can look for correlations of the residuals 
from the relation with other galaxy parameters as a further tool to 
unraveling the processes involved in galaxy evolution. We find that 
there is no correlation of the individual galaxy gradients with any 
other dynamical, morphological or stellar population property dis- 
cussed in this study (excluding the previously mentioned negative 
gradient galaxies), therefore we concern ourselves with examining 
only the offsets from the mean relation. 

While the global residual is an interesting quantity it does not 
take full advantage of the real power of integral-field spectroscopy: 
it's spatially resolved nature. The key result to take from the Mg b- 
V csc relation is that galaxies show the same behaviour both locally 
and globally. The Mg 6-V csc relation holds whether it is measured 
in a central aperture or in annuli. This scale-independent behaviour 
suggests that the Mg 6-V osc relation should hold at all points within 
the central regions of ETGs. It is easy to test this with IFU data by 
plotting the Mg 6-V csc relation for individual bins, completely ig- 
noring any information as to which galaxies any given bin belonged 
to. The result of this analysis is shown in the Fig.QT] Because of the 
sheer number of individual bins in the ATLAS JD data it is neces- 
sary to display this in a different fashion to the elliptical annuli data, 
with the contours reflecting the number of bins at any given point 
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Figure 11. Mgfe-V csc relation determined from individual bins as opposed 
to averaging over elliptical annuli. The solid red line in shows the best-fitting 
relation to the central aperture measurements. The contours of the individual 
bins follow the central relation, with the expected deviation at low V esc due 
to negative gradient galaxies. 

in the Mg 6-V csc plane. As usual the best-fitting global relation is 
overplotted as a solid red line, which follows the contours precisely. 
The bin data are fully consistent with that determined from ellipti- 
cal annuli. The overall behaviour of the relation is also the same, 
with significantly increasing scatter below V csc = 400kms~ . This 
tight relation found even when considering individual bins allows 
us to go further than the globally-averaged mean residual described 
above and consider the residuals of individual bins from the Mg b- 
V csc relation. The individual bin residuals are defined following as 
above but, instead of using global values, using values from indi- 
vidual IFU bins. We will make use of both the global residuals and 
individual IFU bin residuals in the following analysis. 

4.4.1 Deviations from axisymmetry 

As a sanity check on the way our MGE modelling deals with bars 
and other non-axisymmetric features we separate the sample into 
a 'clean' sample of axisymmetric, fast-rotating galaxies that are 
well described by the MGE models and a 'dirty' sample of all 
other galaxies where there are significant deviations between the 
observed photometry and the MGE model. We also include all 
galaxies classified as slow rotators in Paper III as these are likely 
to be weakly triaxial, as well as galaxies with bars, rings, isophotal 
twists and shells or other merger signatures. When we examine the 
Index- V osc relations of these two subsamples separately we find 
that they are consistent with following the same relation, though 
there is larger scatter in the dirty sample. We find o = 0.037 for 
the clean sample, whereas the dirty sample has a = 0.042. The 
increased scatter we find in the dirty sample can account for 15 per 
cent of the scatter in the relation for the complete sample. 

4.4.2 Mg b structures within galaxies 

While we use the same elliptical annuli to measure both the in- 
dex and V csc profiles, the index and V OS c maps often have some- 
what different isocontour shapes. While some of this variation 
can be accounted for by measurement errors there are a num- 
ber of cases where the isocontours of the Mg6 and Fe5015 maps 
show a prominent disk-like structure. This effect is most pro- 
nounced in the Mg b maps. This issue was discussed briefly in 
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Figure 12. Examples of Mg b disk structures in the IFU maps. From top left the galaxies shown are: NGC2577, NGC3377, NGC4473, NGC4570, NGC4638 
and NGC561 1. While the colour scaling is different for each galaxies, the colour bar indicates the general trend in Mg b for all panels The flux contours from 
the collapsed SAURON cubes are overplotted. The Mg b isocontours are flatter than the isophotes in each case. 
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Figure 13. Upper panel: The location of Mg fe-disk galaxies in the Mg b- 
V eS c relation. The Mg b-disk galaxies, indicated by the red profiles, have 
the same distribution as the normal galaxies. Lower panel: Distribution of 
individual galaxy gradients for Mg b-disk (red) and non-disk galaxies (blue). 
The Mg fe-disk galaxies are offset to steeper gradients than non-disk galax- 
ies, however the offset is only marginally significant in a KS-test. 



iKuntschner et alj d2006h andlKrainovic et alj J2008h . and more ex- 
tensively in Kuntschner et al see their fig. 9 and accompany- 
ing discussion). Here we explore how these structures in the index 
maps affect the Mg &-V csc relation. 

We selected by eye all galaxies showing a prominent flattened 
Mgb component in the Mgb maps, yielding a sample of 27 ob- 
jects. A few examples of galaxies with flattened Mgb structures 
are shown in Fig. [T2] These galaxies do not occupy any particu- 
lar region of the Index- V csc diagrams, being indistinguishable from 
early-types that do not exhibit a flattened Mg b component (upper 
panel of Fig. |13t . However, these galaxies have steeper gradients 
than the typical early-type (as well as being steeper than the global 
relation) as shown in the lower panel of Fig.[T3] The shift in slope 
is small, only 1.5<r above the mean local gradient for the whole 
sample. When we exclude the negative gradient galaxies (as these 
are clearly a different population) a Kolmogorov-Smirnov test (KS- 
test) returns the probability that the two sets of gradients are drawn 
from the same distribution is 0.17, which is inconclusive. However, 
these galaxies have residuals from the Mg b-V csc relation consis- 
tent with the full sample of galaxies. This suggests that our choice 
of elliptical apertures is robust against structure in the line strength 
maps. 



4.4.3 Other influences 

In McDermid et al., in preparation we will show that the environ- 
ment and gas content of an early-type galaxy have a small yet sig- 
nificant influence on the global properties of its stellar population. 
These effects persist independent of the size of aperture used to de- 
termine the SSP properties. We briefly revisit these results in this 
paper, using our bin-by-bin analysis and examining the residuals 
from the Mg b-V csc relation (rather than from the SSP-u relations) 
in the context of the local and global behaviour described above. 

In the left panel of Fig. [14] we show the histogram of A Mg b 
for individual bins, divided into bins belonging to all (black), Virgo 
(red) and non- Virgo (blue) galaxies. Bins belonging to Virgo cluster 
galaxies have more positive residuals, in the sense that bins belong- 
ing to Virgo cluster galaxies have higher Mg b at fixed V esc com- 
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Figure 14. Left panel: Histogram of A Mg b for individual bins, divided into all (black), Virgo (red) and non- Virgo (blue) samples. Bins belonging to galaxies 
in the Virgo cluster have higher Mg b at fixed V csc . The Virgo sample does not exhibit the prominent tail of points due to negative gradient galaxies found 
in the full and non- Virgo samples. Right panel: Histogram of A Mg b for individual bins, divided according to the presence of dust (blue) or not (red) in the 
associated galaxy (not in that specific bin). Dusty galaxies tend to have low A Mg b at a given V CS c and also frequently exhibit negative gradients. 



pared to bins belonging to galaxies in lower density environments. 
In the right panel of Fig. [14] we show the histogram of the residu- 
als from the Mg b-V esc of individual bins associated with galaxies 
from the full (black), dusty (blue) and non-dusty (red) samples (as 
classified in Paper II). As can be seen these galaxies are strongly bi- 
ased towards having negative residuals - at a given V eS c a bin from 
a dusty galaxy will typically have a weaker Mg b line-strength than 
a dust-free galaxy. The same result is true of CO and HI detected 
galaxies - bins belonging to galaxies with detectable cold gas com- 
ponents are generally found below the relation. 

These results are consistent with our findings in McDermid et 
al., in preparation. However the use of spatially resolved quantities 
allows us to go one step further in our analysis. As previously stated, 
the Mg b-V C sc gradient does not depend on any galaxy property, in- 
cluding environment and gas and dust content (with the exception 
of the negative gradient galaxies). Here we have shown that individ- 
ual bins within galaxies in high-density environments or with high 
dust or gas contents are offset from the main relation (in a positive 
and negative sense respectively). These two observations show that 
the offset is a systematic one - all bins within the central ~ 1 R e 
are offset from the mean relation by the same amount for a given 
galaxy. This does not apply to the negative gradient galaxies which 
show a broad range of behaviours and no systematic offset. 

4.5 SSP parameters 

In Section POl we explored the connection between the global SSP 
parameters and V csc . In this section we go beyond globally aver- 
aged quantities and instead examine the dependence of the resid- 
uals of the Mg b-V esc relation on th e local SSP parameters. Using 
the SSP models of lSchiavorj J2007l) we apply the SSP method of 
McDermid et al. (in preparation) to the line-strength indices of in- 
dividual IFU bins and construct maps of the three SSP equivalent 
parameters age, metallicity and alpha-enhancement. 

In Fig.[T5]we show the Mg 6-V csc relation for individual bins 
falling within a given range of each of the SSP parameters. The 
coloured contours show the location of all IFU bins falling within 
that range of t (upper panel), [Z/H] (middle panel) or [a/Fe] (lower 
panel) in the Mgb-V eS c plane. The solid line of the same colour 
shows the best fit to the bins with the associated SSP values. The 
contours for each SSP parameter subset are normalised to the total 
number of bins falling within that SSP range - with the contours 



enclosing 10, 30, 50, 70 and 90 per cent of bins with the given SSP 
value. 

For the relation at fixed t or [a/Fe] we find that, except for the 
very youngest or most a-deficient bins, all bins follow the Mg b- 
V csc relation within the 1 a errors. Bins at fixed t or [a/Fe] span 
the full range in V esc of the sample, again except for the youngest 
bins. The youngest bins with t < 1.2 Gyrs (not explicitly shown in 
Fig. 115b all lie below the relation; they are also only found at low 
V csc (consistent with our discussion on negative-gradients galaxies 
in the following Section [5~T| >. Bins with t < 2.4 Gyrs also typically 
lie below the relation, though at higher V CS c this effect is less pro- 
nounced. All bins with t > 2.4 Gyrs follow the Mg b-Y csc relation, 
even at V CS c < 400kms~ 1 . When binning by [a/Fe], only the most 
[a/Fe] deficient bins with [a/Fe] < —0.06 show any significant de- 
viations from the Mg b-Y csc relation (again, not explicitly shown in 
Fig. 1 1 5 1. and only apparent at V CS c < 400kms~ 1 . There is a small 
but significant trend for the relation to be offset to higher values of 
Mg b as [a/Fe] increases: A Mg b increases with increasing [a/Fe]. 
There is no clear trend of A Mg b with t. 

The picture is very different when binning by [Z/H], Bins of 
fixed [Z/H] only span a narrow range in V eS c- The low- [Z/H] bins 
are concentrated at low V esc and lie below the Mg 6-V CS c relation, 
in the region occupied by the negative gradient galaxies. The high- 
[Z/H] bins are concentrated at high V CS c and lie above the main 
Mg b- Vcsc relation. When [Z/H] is restricted to even smaller ranged 
than in the figure, each unique metallicity bin spans only a narrow 
range in V eS c and falls as a clump along the relation. The Mg b- 
Vcsc relation disappears for fixed [Z/H]. This behaviour is not seen 
for t or [a/Fe]; for even narrow ranges of t or [a/Fe] the bins span 
a broad range in V CS c and the Mgb-V CS c relation is still observed. 
This is not surprising, given that the strength of the Mg b absorption 
index is predominantly determined by [Z/H], and so Mg b can vary 
little at fixed [Z/H]. The Mg fe-V eS c relation is predominantly driven 
by the variation of [Z/H] with V eS c, which is again unsurprising 
given the tight mass-metallicity relation of ETGs (McDermid et al, 
in preparation). 



5 DISCUSSION 

In the previous two sections we have presented many pieces of ob- 
servational evidence connected to the role of V OS c in local early- 
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Figure 15. The Mg 6-V csc relation for all individual bins, divided into three ranges of age (upper panel), metallicity (middle panel) and enhancement (lower 
panel). The coloured contours indicate the region occupied by the points with a given range of t, [Z/H] or [a/Fe]. The contours enclose a set fraction of the 
bins with a given range of t, [Z/H] or [a/Fe] - for example there are far fewer total points in the lowest age bin than in intermediate age bins even though the 
same number of contours are displayed. The solid black line shows the global relation for all galaxies while the coloured lines are a fit to the SSP subsets. 
The range of SSP values indicated by each colour is given in the legend on each panel. The negative gradient galaxies stand out clearly at young ages and low 
metallicities. The relation shows the strongest dependence on [Z/H]. 
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type galaxies. In this section we consider the above evidence in con- 
text of galaxy evolution theory and use it to discriminate between 
the different potential formation histories of the nearby early-type 
galaxy population. 

5.1 The effect of recent star formation on the V esc relation 

The most striking aspect of the Mg fe-V esc relation is its tightness 
above V CS c ~ 400kms _1 and the sharp increase in scatter below 
this value. This increase in scatter can largely be attributed to a 
population of galaxies showing negative gradients, with centrally 
depressed values of Mg6 (and to a lesser extent Fe5015) and en- 
hanced values of H/3. These negative gradient galaxies are found 
almost exclusively at low values of V eS c, apart from a single ex- 
ception. As noted above high values of H/3 are usually indicative of 
a young stellar population, with our SSP analysis confirming that 
these galaxies have, on average, lower central SSP ages. All nega- 
tive gradient galaxies have SSP equivalent ages less than 3 Gyrs (5 
Gyrs) within an R e /8 (R e ) aperture (Fig. I16K or, in observational 
terms, H/3 > 2.5(2. 3). In the case of NGC 4150, a prominent nega- 
tive gradient galaxy, ICrockett et alj d201 ll) confirm the presence of 
a centrally concentrated ~ 300 Myr population using HST UV and 
visible photometry. These galaxies also tend to be dusty and contain 
significant amounts of molecular gas, again indicative of recent star 
formation. 75 per cent of negative gradient galaxies are detected 
in CO, compared to 18 per cent of the normal galaxies. Similarly, 
7 1 per cent of negative gradient galaxies contain dust identified in 
the ugriz photometry, compared to only 15 per cent of the normal 
galaxies. 

These young cores typically represent a small (by mass frac- 
tion) perturbation on top of an underlying older population. As the 
galaxy ages and the younger population fades, reducing its contri- 
bution to the total luminosity of the galaxy, we might expect these 
galaxies to return to the Mg 6-V esc relation as the older population 
comes to dominate the light. We can test this idea by artificially 
ageing these negative gradient galaxies using our SSP models. This 
analysis is shown in Fig. [T7] As can be seen the negative gradient 
galaxies typically return to the relation after a few Gyrs of evolu- 
tion, assuming no further star formation takes place. The timescales 
for this are likely shorter than those predicted by the SSP models, 
as the assumptions made in the SSP modelling mean the SSP ages 
are strongly influenced by even a small fraction of very young stars. 
In an upcoming paper in the series we will explore more complex 
star formation histories for the sample, allowing us to accurately 
constrain the mass fractions and ages of these young central stellar 
populations. 

These negative gradients are predominantly found in galax- 
ies with central V osc < 400kms~ 1 . Not all galaxies below this 
Vesc threshold show negative gradients, the majority of galaxies 
in this region still follow the relation. However, these galaxies 
are much more common below this value. Galaxies with central 
Vesc = 400 kms -1 hav e JAM dynamical masses of 3 x 10 10 Mq. 
Kauffmann et al. I d2003l) identified this mass scale as marking the 
transition from young, low surface-density galaxies to old, more 
concentr ated objects. The same value Mjam ~ 3 x 10 1() M s was 
found in lCappellari et al. hereafter Paper XX) to character- 

ize the galaxy distribution on the (M, R c ) and (M, a) planes. This 
mass defines a break in the "zone of exclusion" of the galaxy dis- 
tribution. Below this mass ETGs are found to define a sequence of 
increasing a, bulge fraction, M/L and redder colours. Above this 
mass bulges dominate and the M/L and colour are more homoge- 
neous (see Paper XX for details). 




Age (Gyrs) 



Figure 16. Mg b-W CBC gradient vs. central age, measured within an aperture 
with radius R c /8. The solid line is a fit to the data. The dashed line divides 
galaxies with negative and positive gradients. Negative gradient galaxies all 
have young central ages. 
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Figure 17. The Mg 6-V csc relations for a selection of negative gradient 
galaxies (dashed lines) and their artificially aged counterparts (coloured 
lines). The 'aged' galaxies are produced by keeping their SSP determined 
[Z/H] and [a/Fe] fixed while allowing their SSP age to vary. The colours of 
the 'aged' galaxies indicates the amount of ageing required to return them 
to the Mg 6-V csc relation. Purple/blue corresponds to ~ 2 Gyrs of passive 
ageing, whereas red corresponds to ~ 15 Gyrs. 



5.2 Constraints on merging 

It is perhaps surprising that we do not find any difference, in terms 
of the Mgfc-V C s C relation, between the fast and slow rotators in 
the ATLAS 3D sample, given the differing relative importance of 
gaseous and dissipationless processes in their formation (as dis- 
cussed in Section QJ. Taking an extreme example, if we were to 
merge two completely gas free, equal mass fast rotators that both 
initially lie on the observed Mg &-V csc relation, would we expect 
the resulting merger remnant to lie on the Mg 6-V osc relation? A 
naive initial answer would suggest no; by merging the two galax- 
ies we are significantly altering the mass distribution and hence the 
local V csc , while the stellar population remains unchanged as there 
is no gas to form new stars. If this simple analysis is valid then it 
is possible to constrain the number of major, dry mergers a typi- 
cal early-type galaxy has experienced. However, the picture may be 
considerably more complicated than this simple analysis allows for, 
even in the case of a dry major merger. While the global potential of 
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the galaxy must deepen as more mass is added, the merger process 
may redistribute the stars in such a way that the local V csc is es- 
sentially unchanged. In support of this idea it has long been shown 
that the rank order o f binding energies of the st ars is preserved in a 
collisionless merger dWmtdl978l : lBarneslll988l) . 

We can estimate the effect of a major dry merger by calcu- 
lating the change in V esc from simple analytic expressions. The 
central V csc will be affected by changes in the final mass and ra- 
dius and due to the redistribution of the stellar material due to vio- 
lent re laxation. Following a simple Virial estimate as in lNaab et alj 
J2009h (assuming Virialised initial systems formed dissipatively of 
stars and conservation of energy in the following merger) it can be 
shown that: 

rj_ = (1 + rjf 
n (1 + rje) 

where the subscripts i and / denote initial and final quantities, r 
and M are the radius and mass of the object, r\ is the mass ratio of 
the merger, 1 + rj = Mf /Mi and e is the ratio between the accreted 
and initial mean square speed of the stars, e = (i>a) / (vf).Tf the 
two progenitors are identical (77 = 1, e = 1) the radius and mass of 
the remnant will double. This leads to a change in V esc of: 

Vcsc.f (nMi\ 1/2 = 1 

V csc ,i K \r f Mi J 

For a 1:1 major merger this analytic estimate predicts V OS c will re- 
main unchanged, however this only applies to central values and 
under the assumptions outlined above. If energy is not conserved 
or re laxation is not comp lete the change in V esc may differ. In- 
deed, [Hogkinsleral] I2OO3, their fig. 1) find in their simulations that 

Tj/ri ~ 1.8, as opposed to the value of 2 predicted by the above 

v . f 

analytic estimate. This would yield v ° bc '. ~ 1.05. Equally impor- 
tantly, our analytic estimate is too simplistic to predict how Vesc 
varies with radius after a merger. 

In order to better assess the effect of mergers on the Mg fc-V esc 
relation in real systems we mak e use of the N-body binary merger 
simulations of iBois etal ] J201 j], hereafter Paper VI). By 'tagging' 
the stars in the merger progenitors with a given age and metallicity 
we were able to follow how the stellar population is redistributed in 
a dry major merger. Here we examine only the effect of major (1:1 
mass ratio), dry (< 3 per cent gas) binary mergers on the Mg b- 
V csc relation. In a future work in this series we will fully describe 
this approach including the extension to gas-rich and minor merg- 
ers. We consider two merger cases from Paper VI with early-type 
galaxies as progenitorjfl All models are projected at an inclination 
of 60 degrees to the line of sight. We assign ages and metallicities 
to the progenitor galaxies such that they lie on the observed Mg b- 
V eB c relation. In the model gal axies we use the MILES SSP library 
dSanchez-Blazquez et alj|2006l) to convert the projected metallicity 
and age of the stars into a map of the Mg b line strength. The model 
galaxies are then 'observed' exactly as the real galaxies were and 
Mg b and V osc profiles are constructed as described in Section |4~Tl 

The results of this analysis are shown in Fig.[T8] The progen- 
itors are shown in green, with the corresponding remnant shown 
in blue. The black lines show the observed ATLAS 3D sample and 
the solid red line is the best fitting relation as described in Section 
14.11 The dry merger simulations confirm our naive prediction that 
Mg b remains essentially unchanged while the merger remnants are 

3 their so-called 'remerger 2x11 dd' and 'remerger 2x1 lrr'. See Paper VI 
for details. 
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Figure 18. A zoom in on the high-V eS c end of the Mg b-W CBC relation with 
the results of 1:1 mass-ratio, dry, binary merger simulations overplotted. 
The green lines show the Mg fe-V csc profile for the two (similar) progenitor 
galaxies and the blue line shows the profile for the merger remnant. The 
black lines indicate the observed galaxy profiles. The merger simulations 
are fully described in the text. The effect of the mergers is to increase V eS c 
without significantly changing the observed Mg b. 

shifted to higher V esc by an amount A log V csc = 0.08. We also 
find that Mg £>-V esc gradients are largely preserved in the merger, 
in that the gradient of the remnant is the average of the progeni- 
tors' gradients. It is clear that dry major mergers move galaxies to 
the right of the relation, rather than along the relation. From this re- 
sult we can draw two conclusions. Firstly, massive early-type galax- 
ies are unlikely to have been formed purely by multiple dry major 
mergers of low-mass present-day early-type galaxies. An increase 
in their log V esc by ~ 0.4 due to major mergers (moving them from 
the middle to the high end of the V esc range of our sample), with 
no corresponding change in log Mg b would give these hypotheti- 
cal galaxies values of log Mg b > 3a below the relation. Secondly, 
the observed scatter in the Mg &-V esc relation can be used to con- 
strain the number of major dry mergers a typical early-type galaxy 
can have undergone. Given an observed scatter of a — 0.04, and 
a typical A log V csc of 0.08 per major merger, we find that a typ- 
ical early-type galaxy can only have experienced ~ 1.5 major dry 
mergers since joining the Mg &-V osc relation. 



6 CONCLUSIONS 

In this work we have presented the Index- V csc relations for the 
ATLAS 3D survey, a complete sample of nearby early-type galaxies. 
This survey is based on SAURON integral-field spectroscopy and 
is supported by ugriz photometry. Using these data we accurately 
modelled the mass distributions of 258 galaxies in the sample using 
a combination of MGE photometric models and IAM dynamical 
modelling. We again emphasise that our V csc is derived from the 
surface brightness, so may differ from the true V csc due to the influ- 
ence of dark matter, and that the results we present here are predom- 
inantly empirical rather than causal. An important, and straightfor- 
ward, test of a potential causal relationship would be to examine the 
stellar population - V osc relations are large radii, where dark mat- 
ter is expected to dominate the gravitational potential. From these 
results we are able to draw the following conclusions: 

i) The Mg b - V csc and relation is extremely tight (with a scatter 
of 7 per cent) and displays a local and global connection, in that 
the relation is the same within individual galaxies as it is between 
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different galaxies. The Colour - V csc relations exhibit very similar 
behaviour. These correlations provide a compact way to describe 
gradients in stellar populations in ETGs. They are an ideal tool for 
comparison with numerical simulations, to test whether they repro- 
duce real galaxies. 

ii) Below V csc ~ 400 kms -1 the scatter in the Mg& - V osc rela- 
tion increases significantly due to a population of negative gradient 
galaxies. These negative gradient galaxies typically have centrally 
concentrated young stellar populations and often contain molecular 
gas and dust. This strongly suggests that these negative gradients 
are due to recent star formation. The SSP models suggest that as 
these galaxies age the majority will return to the Mg b - V osc rela- 
tion - negative gradients are largely a transient phenomena. 

iii) Excluding the negative gradient galaxies, the Mg b - V esc gradi- 
ents do not correlate with any other galaxy property. The residuals 
from the V osc relation correlate with: environment, dust and molec- 
ular gas content and a enhancement. 

iv) In the four dimensional space of V esc , Age, [Z/H], [a/Fe] 
the galaxies in our sample are approximately confined to a plane. 
This plane can be defined using just three variables: log V CS c = 
0.62[Z /H] + 0.26 log t + 2.53. This SSP relation is tighter than any 
of the individual SSP - V csc relations and also exhibits the same 
local and global behaviour as the Mg b - V csc relation. 

v) At the high V OS c end of the Mg b - V csc relation a galaxy can 
typically have experienced only ~ 1.5 dry mergers without shifting 
the galaxy so far away from the relation that it is no longer consis- 
tent with the observed scatter. 
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APPENDIX A: MGE MODELS FOR ALL 258 GALAXIES 

Contour plots showing observations (black) and MGE models (red) 
for all 258 galaxies with r-band imaging available. We show only 
the central 40" x 40", approximately the region covered by the 
SAURON observations. The MGE models are able to reproduce a 
range of morphologies from round galaxies to edge-on, dusty, disky 
objects. Bars, which our MGE modelling procedure is designed not 
to fit, stand out clearly. Contours are shown at 0.5 magnitude inter- 
vals down to the standard deviation of the sky background level. 



APPENDIX B: COMPARISON BETWEEN MGE AND 
LITERATURE SURFACE BRIGHTNESS PROFILES 

Comparison of surfa ce brightness profiles d erived from our MGE 
models with those o f iKormendv et al.l d2009h . for all 23 galaxies in 
common between the two samples. The agreement between the two 
sets of profiles is excellent, with typical residuals of 0.04 magni- 
tudes. 
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Figure Al. MGE models (red) and observations (black) for all 258 galaxies with r-band imaging available, arranged by name. 
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Fig. lAll continued. 
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Figure Bl. Comparison of surface brightness profiles from MGE models (solid line) and observations from iKormendv et alj i2009h for all 23 galaxies in 
common between the two samples. The comparison is shown only over the region for which the MGE is valid, 1.5" < r < 2max(<r M Q E ). The sky 
background level in the SDSS r-ban d images is indicated by the horizontal dotted line. 
© 2012 RAS, MNRAS 000. [711191 



